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ABSTRACT 

Extreme mass ratio bursts (EMRBs) have been proposed as a possible source for future space- 
borne gravitational wave detectors, such as the Laser Interferometer Space Antenna (LISA). These 
events are characterized by long-period, nearly-radial orbits of compact objects around a central 
massive black hole. The gravitational radiation emitted during such events consists of a short burst, 
corresponding to periapse passage, followed by a longer, silent interval. In this paper we investigate 
the impact of including relativistic corrections to the description of the compact object's trajectory 
via a geodesic treatment, as well as including higher-order multipole corrections in the waveform 
calculation. The degree to which the relativistic corrections are important depends on the EMRB's 
orbital parameters. We find that relativistic EMRBs (v max /c > 0.25) are not rare and actually account 
for approximately half of the events in our astrophysical model. The relativistic corrections tend to 
significantly change the waveform amplitude and phase relative to a Newtonian description, although 
some of this dephasing could be mimicked by parameter errors. The dephasing over several bursts 
could be of particular importance not only to gravitational wave detection, but also to parameter 
estimation, since it is highly correlated to the spin of the massive black hole. Consequently, we 
postulate that if a relativistic EMRB is detected, such dephasing might be used to probe the relativistic 
character of the massive black hole and obtain information about its spin. 

Subject headings: black hole physics — Galaxy: nucleus — gravitational waves — stellar dynamics 



1. INTRODUCTION 



Low-frequency (1CV 5 < / < 0.1 Hz) gravitational 
wave interferometers, such as the proposed Laser In- 
terferometer Space Antenna (LISA! dBender et al. fl998h 
iDanzmann fcR iidiger 2003; ISumner fc Shau1ll2004f) . will 



open a completely new window to the Universe. Through 
observations of low-frequency gravitational radiation we 
will be able to witness the inspiral and merger of mas- 
sive black hole binaries; the inspiral of compact objects 
into massive black holes; and millions of quasi-stationary 
compact galactic binaries. Recently, a new source of low- 
frequency gravitational radiation has been suggested: ex- 
trem e mass ratio bursts (EMRBs) ([Rubbo et alj 120061 . 
l2007h . 

EMRBs consist of a stellar-mass compact object (SCO) 
orbiting a massive black hole (MBH) of 10 4-8 M Q with 
orbital periods greater than T cut = 3 x 10 4 s. The defin- 
ing orbital period cutoff is derived from LISA's lower 
frequency limit of / cut = 3x 10 -5 Hz. Systems with or- 
bital periods less than T cut will radiate continuously in- 
side the LISA band. Such continuous systems are more 
appropriately categorized as extreme mass ratio inspi- 
rals (EMRIs) and have been studied extensively else- 
where: recent estimations of the event rate are given 
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bv iGair et all (|2004l) and IHopman fc Alexander (|2006T> : 
a discussion on a p o ssible EMRI background is given 
bv iBarack fc Cutlerl (p00 4a): accounts on the theoreti- 
cal desc r iption of E MRIs can be found in the reviews by 
iPoissonl (|2004f ) and iGlampedaklsl (E 005), while the re- 
cent review bv lAmaro-Seoane et al.l ()2007l ) describes the 
astrophysical and detection applications. 

Although EMRB events are distinct from EMRI 
events, their evolutionary track could be connected. In 
the burst scenario, the SCO orbits the MBH emitting a 
beamed burst of gravitational radiation during pericenter 
passage. The emitted radiation carries away energy and 
angular momentum from the system so that after multi- 
ple pericenter passages the orbital period decreases, and 
possibly the system becomes an EMRI. However, this 
evolutionary track is most likely disrupted by scattering 
interactions with other stars and/or if the SCO plunges 
directly into the central MBH on one of its passages. 

The EMRB event rate has recently been investi- 
gated using si mplified galactic models and data analy- 
sis te chniques (jRubbo et alj 120061 120071 : IHopman et al.1 
l2006h. Using a d ensity profile described by an 77- m o del 
(|Tremaine et all 119941 ) . (Rubbo et all (|2006l 120071 ) sug- 
gested an event rate of ~ 15 yr _1 for events with 
signal-to-noise ratios (SNRs) greater than five out to 
the Virgo cluster. When mass segregation and dif- 
ferent inner cusp models are considered, the predicted 
rate decreases by an order of magnitude (jHopman et al.1 
2006). These preliminary studies were aimed at under- 
standing if EMRB event rates are interesting for low- 
frequency gravitational wave detectors such as LISA. 
More work is still needed to improve the predicted event 
rate in the context of realistic galaxies, where the role of 
non-equilibrium dynamics, anisotropy, complex star for- 
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mation histories, substructure, and non-sphericity may 
act to change the rat e from these fiducial estimates by 
orders of magnitude (iHollev-Bockelmann fc Sigurdssonl 
12001 iRubbo et al.ll2006l 120071) . 

In addition to the astrophysical uncertainties, there 
are no investigations of the impact of relativistic cor- 
rections to EMRB dynamics. All EMRB studies have 
been carried out in a quasi- Newtonian framework. In this 
framework, one uses the Newtonian equations of motion 
and extracts the gravitational waveforms by means of 
the quadrupole formula. This approximation ignores the 
black hole nature of the central potential, including the 
black hole's rotation (spin), and is technically valid only 
for orbits with non-relativistic velocities. However, a con- 
siderable number of EMRBs are characterized by large 
pericenter velocities (v p > 0.25 c) and these relativistic 
EMRBs should produce gravitational wave signals with 
larger SNRs, as we will show later. 

In this paper, we shall not consider the issue of EMRB 
event rates, but instead we shall study the effects of rel- 
ativistic corrections to such events. For extreme-mass- 
ratio systems a simple way of introducing relativistic cor- 
rections is by using th e so-called semi-re l ativis tic approx- 
imation introduced bv lRuffini fc Sasaki (11981D, and used 
recen tly in the context of EMRIs by iGair et ail (|2005l 
2006). In this approximation, the MBH and surrounding 
area are modeled using the Kerr solution to Einstein's 
field equations, which describes a stationary spinning 
black hole (the Schwarzschild solution corresponds to the 
non-spinning case). The SCO is considered to be a point- 
like object (neglecting its own self-gravity) whose trajec- 
tory is described by a geodesic of the Kerr spacetime. 
In other words, relative to previous work, we have re- 
placed the Newtonian equations of motion by relativistic 
geodesic equations of motion. 

The relativistic description introduces effects such as 
orbital precession and frame dragging, but it does not 
account for effects due to the gravitational field induced 
by the SCO. These effects, for example, lead to changes 
in the (geodesic) constants of motion due to radiation 
reaction. Even though these effects introdu ce errors that 
scale with the system's mass ratio (e.g. see iGlampedakisI 
2005), they cannot be neglected for EMRIs. This is be- 
cause in the late stages of the EMRI the SCO spends 
a substantial fraction of cycles in the strong-field region 
of the MBH. On the other hand, in the case of EMRBs, 
the SCO sling-shots around the MBH and its interac- 
tion time during pericenter passage is relatively small 
(< 10 5 s). Radiation reaction effects can then be ne- 
glected in EMRBs since the radiation reaction timescale 
is always much larger than the period of pericenter pas- 
sage. 

In this paper, we also improve on the semi-relativistic 
approximation by using a more precise gravitational wave 
extraction procedure. The procedure employed is the 
multipole-moment wave generation formalism for slow- 
motion object s with arbitrarily strong internal gravity 
(|Thornd [l980). We consider terms up to the mass- 
octopole and current-quadrupole multipo les, thus im- 
proving on t he m ass-quadrupole analysi s of IRubbo et al.l 
(|2006L I2007| ) and iHopman et alJ (|2006h . Higher multi- 
poles will become important if the system becomes even 
more relativistic, but pericenter velocities for EMRBs 
are typically small to moderate relative to the speed of 



light (typically 0.1 < v p /c < 0.5). Such higher mul- 
tip olar corrections we re taken into account for EMRIs 
by iBabak et al.l (|2007l ). but for those sources the phase 
evolution must be tracked very accurately, r equiring 
techniques from black ho le perturbation theory (|Poissonl 
12004 lGlampedakidl2005[ ) that we shall not consider here. 

The study of the relativistic corrections considered in 
this work leads to the following conclusions. First, we 
find that relativistic effects are significant for approxi- 
mately 50% of the o rbits containe d in th e EMRB phase 
space considered by IRubbo et all (|2006l I2007D . These 
relativistic EMRB orbits differ from their Newtonian 
counterparts in such a way that the associated waveforms 
present a noticeably different structure. In particular, 
we find that there is a dephasing relative to Newtonian 
waveforms that is due to precessional effects and depend 
strongly on the MBH spin. These findings show that 
EMRB events are relativistic enough that they should 
be t reated accordingly, as was previously found for EM- 
RIs (IGlampedakisI l2005h . 

Second, we find that the corrections to the trajecto- 
ries affect the waveforms much more than the correc- 
tions in the waveform generation over several bursts. 
For example, for a given relativistic trajectory, we find 
that the difference between the SNR of a waveform ob- 
tain from the quadrupole formula to that obtained from 
the quadrupole-octopole formula is of the order of 10% 
(depending on the location of the observer.) On the 
other hand, using the same waveform generation formula 
(quadrupole or quadrupole-octopole), the difference be- 
tween the SNR of a Newtonian waveform to that of a 
Kerr waveform is of the order of 100%. These find- 
ings show that modeling EMRB waveforms with a quasi- 
Newtonian treatment might not be sufficient for certain 
data analysis applications. 

Third, we find that the relativistic corrections accu- 
mulate with multiple bursts and, thus, they may have 
an important impact in improving the SNR. It is also 
conceivable that such corrections might be important for 
parameter estimation studies and, perhaps, may be used 
to determine or bound the spin of the MBH if a high 
SNR event is detected. Along this same lines, if parame- 
ters can be determined accurately enough, it might also 
be possible to use EMRB measurements to test devia- 
tions from General Relativity. We must note, however, 
that changing the orbital parameters in Newtonian wave- 
forms could somewhat mimic some of the relativistic cor- 
rections, but a detailed Fisher analysis of such effects is 
beyond the scope of this paper. 

The remainder of this paper is divided as follows: Sec- 
tion [2] deals with the dynamics of EMRBs in the semi- 
relativistic approximation and justifies the use of this 
approximation for these systems; Section [3] reviews the 
inclusion of higher-order multipolar corrections to the 
waveform generation formalism; Section |4] describes the 
numerical implementation of the equations of motion and 
the initial data used; Section [5] compares the orbital tra- 
jectories and waveforms; Section [6] concludes and points 
to future research. 

In this paper, we denote the MBH mass by M, and 
its qravitational radius by R, — 2GM,/c 2 , where c is the 
speed of light and G the Newtonian gravitational con- 
stant. To simplify some expressions we normalize masses 
with respect to M MW = 4 x 10 6 M , which is of the same 
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magnitude as the mass of the MBH at the center of the 
Milky Way (|Ghez et al.ll2005h . The gravitational radius 
can then be written as: 



R. = (3.82 x 1CT 7 pc) 



M. 



(1) 



2. EMRB DYNAMICS 



In this section, we discuss the description of the orbital 
motion. Newtonian dynamics usually provides an ade- 
quate description of many astrophysical sources of grav- 
itational waves, at least from a qualitative point of view. 
However, for certain gravitational wave sources, such a 
description is insufficient and relativistic effects have to 
be considered. For EMRB sources with pericenter dis- 
tances r p > 4i?. and velocities v p /c < 0.5, the semi- 
relativistic approximation to the equations of motion, in 
combination with a multipolar description of the grav- 
itational radiation, can adequately model the dynamics 
and gravitational radiation, as we argue below. 

The semi-relativistic approximation treats the motion 
of the SCO in the point-particle limit as a geodesic of 
the MBH geometry, which is justified based on the small 
mass-ratios associated with these systems. In this work, 
we adopt Cartesian Kerr-Schild coordinates, {t, x 1 } (i = 
1, 2, 3), in which the MBH geometry is time-independent, 
reflecting its stationary character, and tends to a flat- 
space geometry in Cartesian coordinates far from the 
MBH. We denote the geodesic trajectory by z l (t), its 
spatial velocity by u l (t) = dz % /dt, and its spatial acceler- 
ation by a l {t) — dv l /dt. The latter, in such a coordinate 
system, and by virtue of the geodes ic equations of mo- 
tion, has the following form (e.g. see lMarck1 [l996): 



(2) 



where g„ y (p,v = 0,1,2,3) are the spacetime compo- 
nents of the MBH metric. These equations describe the 
influence of the spacetime curvature produced by the 
MBH and approach the Newtonian equations of motion 
in the regime where v/c= \v l \/c <C 1 and GM,/(c 2 r) <C 
1 (r=\x% 

The effects from the self-gravity of the SCO can be 
neglected. To see this, consider the (Keplerian) or- 
bital timescale, T orb , in comparison to a characteristic 
radiation-reaction timescale, T rr . For the latter, we can 
use the timescale associated with the rate of change of the 
semi-latus rectum, p, related to the pericenter distance 
byr p = p/(l+e), namely T rr ~ p/\dp/dt\. The radiation- 
reaction timescales of the other orbital elemen ts are com- 
parable or larger (see, e.g.. lGlampedak is 2005). The ratio 
of these timescales is 



T. 



2nfi 



m 



5/2 



(3) 



where /i = m /M m is the mass ratio of the system and m 
the SCO mass. It is evident that the radiation-reaction 
timescale is much greater than the orbital timescale due 
to the extreme mass ratio, fi <C 1, and because EM- 
RBs have p ^> fi 5 / 2 R,. In the unlikely case that more 
accuracy is required, one could improve the analysis 
throu gh the use of "Kludge" waveforms (iBabak et al.l 
I2007T) , which have been shown to reproduce numerical re- 
sults in the adiabatic approximation accurately for EM- 
RIs. 



Formally, the orbital timescale used is not really the 
exact timescale of orbital motion. This is because the 
mass distribution of a MBH-embedded galaxy possesses 
a non-Keplerian potential that leads to non-Keplerian 
orbits. However, most EMRBs (by rate) have apocen- 
ters that do not extend far into the stellar population, 
implying that the contribution from the galaxy potential 
is minimal. The orbits we study in later sections have 
a contamination from the galactic poten tial that is less 
than 2% of the MBH mass. Moreover, Hopman et 
(2006) rightly argue that the inner region is statistically 
empty of stars, which is due to finite effects realized at 
the small scales observed near the MBH. 

Certain constraints may be derived on the size of p and 
r p for EMRB events. The most important constraint 
is derived from the definition of EMRBs: orbits with 
sufficiently large orbital period T or b > T cut . Assuming 
a Keplerian orbit (which is a rough assumption), this 
constraint translates to pericenter distances as follows 



r p > (7.98 x 10~ 7 pc) 



(1-e) / M. 
0.1 \M Mn 



1/3 



T 

T C ut 



2/3 



(4) 

where we have resca led quantities assuming a typical ec- 
centricity of e = 0.9 (|Rubbo et al.ll2006Ll2007h and a typ- 
ical MBH mass of M, = M MW . In terms of geometrized 
units, such a constraint translates roughly to r p > 2R,. 

This constraint can be compared with the requirement 
that the SCO does not get captured. Any object that 
enters the black hole event horizon is captured, where 
the horizon is located (in Boyer-Lindquist coordinates) 
at 



(1.91 x 10" 7 pc) 



M. 




, (5) 



where a. is the (Kerr) MBH spin parameter, related to 
its intrinsic angular momentum by S, — GM, a,/c, and 
bounded by a, /ill, < 1. Thus, for a maximally spin- 
ning Kerr MBH (a. = A/.), r cap = 0.5i?., while for a 
Schwarzschild (non-spinning) MBH it is just R,. This 
condition tell us simply that r p > r cap , which is a condi- 
tion superseded by the constraint on the orbital period 
given in equati on (|4l). One could expl ore other possi- 
ble constraints (|Rubbo et alJl2006l 12007) but they are in 
general superseded by equation (TJJ. 

These constraints clearly exclude the ergosphere of the 
MBH (r cap < r < R,) where frame dragging effects are 
most pronounced. Howev er, for EMRIs it has been ar- 
gued (Glampedakis] l2005f ) that orbits with r p < 10R, 
cannot be considered Keplerian anymore, mainly due to 
precessional effects. This statement can be made more 
quantitative by looking at the ratio of first-or der post- 
Newtonian (1 PN) predictions (|Blanchetll2006T ) to New- 
tonian ones (0 PN). For example, for the energy of a cir- 
cular orbit, this ratio scales as 7i?./(8r p ), while for the 
perihelion precession angle, the ratio scales as 3i?./ (2r p ), 
for extreme-mass ratios. Therefore, for orbits with peri- 
center passage r p ~ 5i?., the 1 PN correction to the 
energy and the perihelion precession angle is approxi- 
mately 20% and 30% respectively, relative to the New- 
tonian value. This indicates that, even for orbits outside 
the ergosphere, relativistic effects are not necessarily neg- 
ligible. 
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The relativistic geodesic equations of motion introduce 
corrections to the Newtonian motion that can be inter- 
preted in terms of a black hole effective potential. By 
comparing the Newtonian and relativistic potentials one 
can see that the relativistic corrections dominate over 
the centrifugal barrier at small distances from the black 
hole center. In this work we show that these relativistic 
corrections can be sampled by EMRBs and hence, one 
should model these systems accordingly. Nevertheless, 
as we argued above, the relativistic treatment of EMRBs 
does not need to be as sophisticated as in the case of 
EMRIs, since radiation-reaction can be neglected. 

3. EMRB WAVEFORMS 

In this section we describe how we extract gravita- 
tional waveforms once we have integrated the geodesic 
equations of motion. We use a multipole-moment 
wave generation formalism for slowly- moving objects 
with arbitrarily strong internal gravity (iThornd 1 19801 ; 
iFlanagan fc Hughesll2005t lGlampedakiil2005l ) . In quasi- 
Newtonian and semi-relativistic treatments, the radia- 
tion is modeled by the lowest non-vanishing multipole 
moment: the mass-quadrupole. To that order, and for 
the case of a point-like object orbiting a MBH at a fixed 
coordinat e location, the plus and cross po larizations are 
given by (jMisner et al.lll973t lThornd[l980l ) 



_ 2Gm i:j 

H>+,X — 4 e +,x \ a i z j 



(6) 



where r is the (flat-space) distance to the observer and 
el_ x are polarization tensors. This expression assumes, 
based on the slow motion approximation, that the change 
in the acceleration with respect to time, the jerk, j l = 
da 1 /dt, is a small quantity. More precisely, we are ne- 
glecting terms of order (v/c) 3 , or in other words, since 
the (quadrupole) leading order terms are of order (v/c) 2 , 
this implies a relative error of order v/c. 

One can improve on this description for the gravita- 
tional radiation by accounting for higher-order multipole 
moments. In this paper, we consider the mass-octopole 
and current-quadrupolc multipoles, which require the 
knowledge of one more time derivative of the trajectory 
the jerk. Adding these contributions t o equation Q , the 
gravitational waveforms are given by (|Thorneill980f) 

2Gm f j 



— C+)X >, ciiZj + ViVj 

(n • z) (zijj + SaiVj) + (n • v) (ojZj + v^j) 
1 



(n • a) ViZj - - (n ■ j) ZiZj 



(7) 



where n l = x l /r is a unit vector that points to the ob- 
server and the vector product is the flat-space scalar 
product. In this case, we are neglecting terms of order 
{v/c) 4 and hence we are making a relative error of or- 
der (v/c) 2 with respect to the leading order quadrupole 
term. 

The waveforms of equation |(7J) are a truncated mul- 
tipole expansion, where we are neglecting the current- 
octopole, mass-hexadecapole, and higher multipole mo- 
ments. This expansion is based on a slow-motion approx- 
imation which is valid for orbits whose pericenter velocity 



is small relative to the speed of light. For closed circu- 
lar orbits, we can use the Virial theorem to argue that 
this is equivalent to requiring r p > M m . For a relativistic 
EMRB event with v p /c — 0.4, the maximum relative con- 
tribution of the octopole to the quadrupole is of the order 
of 40%, since the octopolar term is of order v/c smaller 
than the quadrupolar one. In this paper, we shall study 
E MRBs from the sam ple of Milky Way sources studied 
in iRubbo et al.l (|2006f ) . These sources have initial peri- 
center distances of r p > 8M m , thus justifying the use of a 
low-multipolar expansion in the wave-generation formal- 
ism and the neglect of radiation reaction effects in the 
orbital motion. 

4. NUMERICAL SIMULATIONS 

In this section we describe the EMRB simulations that 
were carried out, including the choice of initial condi- 
tions. The simulations involve integrating the equations 
for geodesic motion around a Kerr black hole, equa- 
tion @, forward in time. (For a deta iled exposition of 
Kerr geodesies see lChandrasekharll992D . Since Cartesian 
Kerr-Schild coordinates are used, the initial conditions 
can be denoted by (zq,Vq). The numerical implementa- 
tion does not use the Kerr geodesic constants of motion 
(energy, angular momentum, and Carter constant) in or- 
der to reduce the number of variables of the resulting 
system of ordinary differential equations. Instead, we 
have used the constants of motion to monitor the accu- 
racy of the time integration. The integration accuracy is 
set so that we obtain fractional errors for the constants 
of motion smaller than one part in 10 10 . The code uses 
a Bulirsh-Stoer extr apolation method as the evolution 
algor ithm (see, e.g. iPress et ail Il992t iStoer fc Bulirschl 
1993). We have also introduced in the code the possibil- 
ity of switching between Kerr geodesies and Newtonian 
equations of motion. The gravitational waveforms are 
then obtained directly by applying expressions © and 
to the numerically obtained trajectory z l (t). 

Comparisons are carried out by choosing a represen- 
tative relativistic orbit within the allowed phase space 
for EMRB's. We made the following choices for the test 
case: 

• The central MBH mass is M. = M MW and the SCO 

mass is rn = 1 Mq, such that the mass ratio is 
fi = m/M. = 2.5 x 10" 7 < 1. 

• The MBH spin parameter is either a. = 
(Schwarzschild) or a. = 0.998 M, (Kerr). The an- 
gular momentum is aligned along the z-axis and 
equal to either S z = or S z = 0.998GM 2 /c. 

• The observer is located at r (, s = 8 kpc along the 
z-axis, which corresponds to the approximate dis- 
tance from Earth to th e center of the Milky Way 
(|Eisenhauer et al1l2005f) . 

Furthermore, we make the following choices for the or- 
bital initial conditions: 

z l Q = (-1.59, 1.05, -0.185) x 10~ 5 pc, 

v l = (1.70, -2.89, 0.510) x 10 4 km s _1 . (8) 

The initial conditions are such that tq = \zq\ = 50i? # = 
1.91 x 10~ 5 pc, and \v* \ = 0.11 c = 3.39 x 10 4 km s _1 . 
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The orbital plane is inclined by 10° with respect to the 
x — y plane to demonstrate the effects of orbital plane 
precession, which only occurs for spinning MBHs. Since 
this paper is concerned with the effect of relativistic cor- 
rections to EMRB events, we choose to give all orbits the 
same initial conditions. The possibility of relaxing this 
choice and its effect on the conclusions derived in this 
paper shall be discussed in a later section. 

Although the t est orbit is in the pha se space of EMRB 
events studied in IRubbo et al.l (|2006f) . one might worry 
that it is too relativistic to actually have a significant 
probability to occur in nature. In particular, one can 
think that the SCO may be tidally disrupted. To address 
this question, let us consider a Newtonian description of 
the central potential, which leads to the following values 
of the pericenter distance and velocity: 

r p = 6.45 R. = 2.48 x 1CT 6 pc, 
| Up | = 0.384 c= 1.15 x 10 5 km s" 1 . (9) 

One might worry that stars might be tidally disrupted 
with such sma l l per icenters. However, as shown by 
iHopman et~aH (|2006D . most SCOs in EMRB scenar- 
ios consist of stellar-mass black holes, which cannot be 
tidally disrupted. 

Such relativistic orbits are actually naturally occur- 
rin g in the ph a se spa ce of possible EMRBs studied 
in IRubbo et all (|2006h . Of all o rbits in the a l lowed 
EMR B phase space considered in IRubbo et al.l (2006, 
l2007f) . 6% is contained within a small 6-dimensional 
phase space volume centered on the test orbit 6 . Fur- 
thermore, the test case shown possesses a short orbital 
timescale, bursting approximately 150 times per year. 
Such events with small orbital time scale were shown t o 
dominate the EMRB event rate in IRubbo et al.l (|2006f ). 
It is in this sense that the test orbit studied here is typical 
or representative of EMRBs. 

The relative location of the test orbit in the pericen- 
ter distance-ecc entricity plane of t he phase space of al- 
lowed EMRBs (jRubbo et al.ll2006l ) is presented in Fig- 
ure [T] (triangle.) The eccentricity was here calculated as- 
suming a Newtonian orbit and the pericenter separation 
is given in gravitational radii, i?,. Although the test or- 
bit has a large eccentricity, its apocenter is small enough 
i r a 150i?. w 6 x 10 -5 pc) that the contribution from 
the surrounding stellar population to the potential can 
be neglected. In general, the left side of the figure corre- 
sponds to highly relativistic orbits with large pericenter 
velocities and small pericenter distances. Orbits with 
pericenter velocities \v p \ > 0.25 c = 3 x 10 4 km s _1 ac- 
count for approximately 50 % of the possible orb its within 
the phase space studied in IRubbo et alJ (|2006[ L 

Although the test orbit is a good source to demon- 
strate the differences between the Newtonian and rela- 
tivistic treatments, we could have chosen an even more 
relativistic one. Such an event would still be classified 
as an EMRB in a Newtonian treatment, but it would 
border with the definition of a continuous source. An 
example of such an extreme orbit is shown with a circle 

6 In other words, 6% of the test orbits considered to be EMRBs in 
IRubbo et al. ( 200g) are close in phase space to our test orbit. This 
does not imply that the probability of such a test orbit actually 
occurring in nature is 6%, since any element of phase space may 
have small overall probabilities, even down to 10 _7 %. 




Fig. 1. — Plot of th e set of possible EMRB orbits as computed 
in lRubbo et al.l l|2006) in the pericenter distance-eccentricity plane. 
The pericenter distance is given in units of the gravitational radius 
R, . The initial conditions for the test [Eq. {8}] and the extreme 
[Eq. I|10|l ] orbits are indicated by a triangle and a circle respectively. 

in Figure [U to the left of the test orbit (triangle). This 
extreme orbit possesses the following initial conditions: 

4 = (-1.81, 0.6, -1.06) x 10~ 6 pc, 

v' = (1.72, -1.78, 0.31) x 10 5 km s"\ (10) 

where r p = 4R. = 1.53 x 10~ 6 pc and \v p \ = 0.49c = 
1.46 x 10 5 km s _1 for a Newtonian potential. We will 
study such an extreme orbit at the end of the next section 
as an example of a limiting relativistic case. 

5. COMPARISON OF TRAJECTORIES AND WAVEFORMS 

In this section we compare the results obtained for 
both the orbital motion and the gravitational radiation 
emitted by an EMRB event using both the Newtonian 
and relativistic description. Since the plus and cross po- 
larization waveforms present similar features, we only 
plot the plus polarization waveforms. In the remain- 
der of this section we use the following nomenclature: 
a quadrupolar (octopolar) Newtonian waveform is one 
that was calculated using the quadrupole (octopole) for- 
mula and Newtonian equations of motion; a quadrupole 
(octopole) Schwarzschild waveform is one that was cal- 
culated using the quadrupole (octopole) formula and the 
geodesic equations of motion with no spin (a. = 0); a 
quadrupole (octopole) Kerr waveform is one that was cal- 
culated using the quadrupole (octopole) formula and the 
geodesic equations of motion with spin a. = 0.998M.. 

5.1. Orbital Trajectories 

Let us begin by comparing the trajectories obtained 
in our simulations. In Figure [2 we plot the test orbit 
corresponding to a Newtonian treatment (solid line) and 
the one corresponding to a relativistic treatment with- 
out spin (dashed line) and with spin (dotted line). The 
dot and arrow indicate the initial location and velocity 
projected onto the x-y plane. The MBH is located at 
the origin of the coordinates, and the vectors denoted 
by L and S describe the direction of the initial orbital 
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Fig. 2. — Trajectories for the SCO, with initial conditions given 
by J5J, corresponding to a Newtonian description (solid line) and 
rclativistic descriptions with no spin (dashed line) and with spin 
(dotted line). The MBH is located at the origin and the vectors L 
and S denote the initial orbital angular momentum and the MBH 
spin respectively. 

angular momentum and the MBH spin respectively. In 
the relativistic description there are precessional effects 
in the SCO orbit that can be clearly observed in Fig- 
ure O These precessional effects are: pericenter preces- 
sion about the orbital angular momentum axis, which 
acts in the initial orbital plane; and frame-dragging pre- 
cession about the total angular momentum axis, which 
acts out of the initial orbital plane. While the former 
always occurs in a relativistic treatment, the latter is 
present only in the spinning case. 



Different relativistic precessional effects are generally 
of different magnitude. These effects are usually inversely 
proportional to the pericenter distance, or equivalently 
proportional to the pericenter velocity of the SCO. Pre- 
cession out of the initial orbital plane, however, is smaller 
than precession in the orbital plane by a relative factor 
of order v p /c and it is directly proportional to the spin of 
the MBH. In t erms of post-Newtonian theory (see, e.g. 
lBlanchetl l2006). the pericenter advance is described by 
lst-order post-Newtonian corrections to the equations of 
motion (order {v/c) 2 relative to the Newtonian acceler- 
ation), while precession off the orbital plane is due to 
spin-orbit and spin-spin couplings that correspond to 1.5 
and 2-order corrections (order (v/c) 3 and {v/c) 4 relative 
to the Newtonian acceleration.) Therefore, since EM- 
RBs are characterized as events with small to moderate 
pericenter velocities, precession out of the initial orbital 
plane is small to moderate relative to pericenter advance, 
even for maximally spinning MBHs. 

We can estimate the precession rate by comparing the 
Newtonian and relativistic trajectories. For the test or- 
bit considered, we find that the rate in the orbital plane 
is roughly 7r/3 radians per orbit for the non-spinning 
case and 2ir/3 radians per orbit for the spinning one. 
These precessional effects have b een studied ex tensively 
in the context of EMRIs (see, e.g. Schmidt 2002) and also 
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Fig. 3. — EMRB waveforms (plus polarization): The Newtonian 
waveform corresponds to the solid line and the Schwarzschild one 
to the dashed line. 

sp ecifically for S-st ars in the central region of our Galaxy 
in iKraniot is (2007). Nonetheless, these effects have not 
been previously analyzed in the context of EMRBs, since 
previous studies employed a quasi-Newtonian treatment. 

5.2. Waveforms 

Let us now analyze how the differences in the SCO tra- 
jectories translate into different signatures in the wave- 
forms. In Figure [3] we plot the quadrupole Newtonian 
and Schwarzschild waveforms (solid and dashed lines re- 
spectively), while in Figure 2] we plot the quadrupole 
Schwarzschild and Kerr waveforms (dashed and dot-dot- 
dashed lines respectively.) There are three main differ- 
ences between the Newtonian and the relativistic wave- 
forms: a modulated phase change, an amplitude change, 
and a time of arrival change. The changes in amplitude 
and time of arrival are due to the test particle experi- 
encing a larger "force" of attraction as it approaches the 
black hole. Quantitatively, this increase in force is due 
to the presence of (r~ ")-contributions to the relativistic 
corrections to the central potential (with n a real positive 
number.) 

Gravitational wave interferometers are most sensitive 
to the phase, which is clearly different for Newtonian 
and relativistic waveforms. The dephasing present in 
Figures [3] and Q] parallels the orbital dephasing discussed 
earlier (the gravitational- wave and orbital frequencies are 
intimately related) and leads to an amplitude modula- 
tion. In terms of the gravitational wave phase, both 
Figures [3] and H] show a dephasing of 7r/6 radians per 
cycle. This can be seen after the third burst where the 
waveforms are back in phase. In fact, there is a signif- 
icant dephasing even between the relativistic waveforms 
due to the effect of the MBH spin. If a cursory exami- 
nation by eye can detect the difference in the waveforms 
due to differences in the nature of the central potential, 
it stands to reason that strong-field EMRB waveforms 
might allow us to probe of the spacetime near a MBH. 
We should note, however, that no work has yet been done 
to find best-fit parameters for Newtonian waveforms that 
maximize the correlation with relativistic ones. In other 
words, it might be possible to mimic some of the rela- 
tivistic corrections by choosing different initial data for 
the Newtonian waveforms, but such a study is beyond 
the scope of this paper. 
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Fig. 4. — EMRB waveforms (plus polarization): The 
Schwarzschild waveform corresponds to the dashed line and the 
Kerr one to the dot-dot-dashed line. The dephasing of the wave- 
forms can be best observed during the silent transitions between 
bursts. For example, in the first silent transition the waveforms 
are roughly n radians out of phase, while in the third one they are 
in phase. 

The difference in dephasing can be better studied by 
calculating the signal overlap, 



(hx\h 2 ) = 



fo hi(t)hn(t)dt 



yj 'fi h\(t)dt yj / T h\(t)dt 



(11) 



The overlap indicates how well a signal h\ can be ex- 
tracted via matched filtering with a template h 2 7 ■ In 
Figures [5] and [5] we plot the normalized integrand as a 
function of time, with h\ given by the quadrupole Kerr 
waveform and h 2 given by either the quadrupole Newto- 
nian or Schwarzschild waveforms. Observe that neither 
the Newtonian nor the Schwarzschild waveforms match 
well with the Kerr waveform. Moreover, note that the 
correlation with the Newtonian waveform deteriorates 
greatly after only the first cycle. The integral of equa- 
tion (fTTj) gives the correlation between different wave- 
forms over nine days (four bursts): for the Newtonian 
and Kerr plus-polarized waveforms it is 9.6%; for the 
Schwarzschild and Kerr plus-polarized waveforms it is 
—6.3%. As a point of comparison, a substantial signal 
overlap should be > 95%. As already mentioned, the 
same initial conditions were chosen for both the New- 
tonian and relativistic orbits, such that their waveforms 
would be both initially in-phase and any dephasing due 
to relativistic effects could be clearly studied. However, 
such a choice forces the SCO to pass through periap- 
sis at slightly different times, because in the relativistic 
case this object experiences a "deeper" potential. Such a 
difference in timing degrades the overlap somewhat and 
could, in principle, be ameliorated by choosing different 
initial conditions for the Newtonian evolution, i. e. by 
maximizing the overlap over all orbital parameters, but 
such a study is beyond the scope of this paper. 

Figures [5] and [5] provide some evidence that the use of 
a relativistic waveform might be required for the data 
analysis problem of extracting EMRB signals. Such ex- 
pectations are somewhat confirmed in Table [1] where 

7 The overlap is given here in the time domain, but an anal- 
ogous representation in the frequency domain could also be used. 
Such an expression in the frequency domain can be derived through 
Parseval's theorem. 
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Fig. 5. — Plot of the overlap integrand of equation JTTJ with hi 
given by the quadrupole Kerr waveform and Y12 by the quadrupole 
Newtonian one. The inset zooms to a region near the small peaks. 

we present the correlation between Newtonian and Kerr 
plus-polarized waveforms integrated over both a single 
day (one burst) and nine days of data (several bursts) for 
a sample of different EMRB orbits 8 . Both a single and 
several bursts should be studied because, in principle, 
parameter adjustments might mitigate the between-burst 
dephasing, but probably not the in-burst dephasing. The 
orbits chosen were taken directly fr om the allowed EMRB 
phase space of iRubbo et al.l (|2006h and possess different 
initial positions and velocities, leading to different eccen- 
tricities, initial inclination and orbital periods. The or- 
bital period can be used to classify the orbits into highly- 
bursting (burst more than once per week) and slowly- 
bursting (burst less than once per week.) For highly- 
bursting EMRBs, the average correlation is ~ 27% when 
all nine days of data (several bursts) are used, while 
it is ~ 93% when only one day of data (single burst) 
is used. For the slowly-bursting EMRBs we considered 
here, there is actually only one burst per week and its 
correlation is ~ 85%. These results indicate that accu- 
mulating precession effects lead to a significant loss of 
correlation between Newtonian and Kerr waveforms if 
bursts are to be connected. Moreover, we see that even 
for a single burst the shape of the relativistic waveforms 
is sufficiently different from its Newtonian counterpart to 
lead to a significant mismatch (in-phase dephasing). If 
one maximizes the correlation over intrinsic orbital pa- 
rameters it might be possible to increase the correlation 
somewhat, but again that is to be studied elsewhere. 

Let us now focus on the differences in the waveforms 
when they are calculated with the quadrupolar approx- 
imation versus the quadrupolar-octopolar one. In Fig- 
ure [7] we plot the absolute value of the difference between 
the octopole and quadrupole waveforms as a function of 
time for a Schwarzschild (upper panel) and a Kerr (lower 
panel) central potential. The difference is normalized to 
the maximum of the first peak of the quadrupole wave- 
form, since other peaks have approximately the same 
maximum. Observe that the inclusion of higher-order 
multipoles does not affect the phasing of the waveforms, 
but only the amplitude, which is in general different by 
roughly 4% relative to the quadrupole waveform for the 
Kerr test case. At first sight, this result is in disagree- 

8 Here we use the frequency representation of the correlation 
calculation, employing the Fourier transform of the waveforms. 
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Fig. 7. — Plot of the absolute magnitude of the difference be- 
tween the quadrupole and octopole Schwarzschild (upper panel) 
and Kerr (lower panel) waveforms. The difference is normalized to 
the maximum value of the first peak of the quadrupole waveform. 

ment with the expectation that the octopolar correction 
is at most < 40% of the quadrupolar one. Note, how- 
ever, that the 40% estimate is an order-of-magnitude up- 
per limit, since the octopole correction is dependent on 
the location of the observer relative to the trajectory, ve- 
locity, acceleration and jerk vectors. For the test case, 
where the observer is located on the z-axis and the or- 
bit is initially inclined by 10°, the octopolar change is 
reduced by approximately an order of magnitude, since 
initially (n • z) m (n ■ v) w (n ■ a) fa (n ■ j) « 0.1. In 
Table [2] we present the approximate maximum differ- 
ence between octopolar and quadrupolar waveforms as a 
function of observer location, focusing only on the first 
burst of radiation. The location of the observer is rotated 
about the y-axis on the x — z plane, always at a fixed ra- 
dial distance of 8 kpc. Note that for some observers the 
difference is larger and, in fact, of the order of 40%, since 
the dot products are closer to unity. These results are 
thus consistent with the expectation that the n-th multi- 
polar contribution cannot in general be larger than order 
(v/c) n relative to the quadrupolar leading term. 

5.3. Data Analysis 

In order to quantify some of our statements about the 
change in phase and amplitude, we calculated the SNR 



= 4 



IM/)I 2 

Sn(f) 



df 



(12) 



where the tilde denotes the Fourier transform and 
S n (f) is the one-sided power spectral noise density. 
Her e we employ the O nline Sensitivity Curve Genera- 
tor ([Larson et al. 2000) with the standard LISA settings 
and the inclusion of the white-dwarf background contri- 
bution. When calculating SNRs, we set the observation 
time to roughly nine days, so as to include multiple bursts 
in our single SNR value. 

The inclusion of relativistic corrections in the trajecto- 
ries has a dramatic impact in the SNR. We find that the 
Schwarzschild waveform increases the SNR by a factor of 
approximately 59%, while the Kerr waveform increases 
it by 162%, relative to the Newtonian SNR. The differ- 
ence in SNR is because the relativistic orbits experience 
a deeper effective potential and, thus, the interaction 
timcscale is smaller. Therefore, the inverse of the inter- 
action time, /* = v p /r p , is larger for the Schwarzschild 
and Kerr waveforms relative to the Newtonian one. As a 
result, the Fourier power is shifted to higher frequencies, 
where LISA is more sensitive. 

The SNR behaves similarly for other EMRB orbits 
with different orbital periods, eccentricities and percen- 
ter parameters. This can be observed in Table [T] where 
we present the SNR difference between Newtonian and 
Kerr waveforms for different EMRB orbits for a single 
and several bursts. Highly-bursting EMRBs lead to a 
large change in the SNR over a week of data, since they 
experience the depths of the effective potential several 
times. Per burst, the change in SNR can range from 
one, to ten or even ninety percent, depending on the in- 
clination angle of the orbit, the pericenter distance and 
other orbital parameters. Also note that the orbits pre- 
sented in the table are not as relativistic as the test case, 
which is why the SNR difference is smaller. This study 
seems to indicate that the SNR is in general somewhat 
larger for relativistic waveforms, specially when several 
burst are taken int o account. 9 Consequently , the event 
rate calculated in iRubbo et al.l (|2006l . l2007t ). which in 
particular summed over all bursts in one year of data, is 
an underestimate for their galaxy model, because some 
of the systems with a Newtonian SNR < 5 should have 
been added to the detectable event rate. However, the 
uncertainty in the event rate is still probably dominated 
by astrophysical uncertainties and not by the dynamics 
modeling. 

Conversely, the inclusion of higher multipole moments 
to the wave generation formalism has little to no ef- 
fect in the SNR. In the previous section we showed 
that there was ~ 4% difference between the octopole 
and quadrupole waveforms relative to the maximum of 
the first peak of the quadrupole waveform. We further 
showed that this difference depends on the location of 

9 Naive intuition might suggest that the change in the SNR 
should scale like the square root of the number of bursts, but 
this is not necessarily correct. First, different orbits possess dif- 
ferent beaming patterns on the sky due to precession. Second, 
the starting frequencies of these bursts is very close to the limit 
of LISA'S sensitivity band [10 Hz], As the orbits burst, pre- 
cession somewhat increases the frequency of the waveform, forcing 
different bursts to contribute different amounts to the SNR. 
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the observer (see table [2]), but for the test case it does 
not exceed a maximum of 40%, which agrees with the 
multipole-ordering argument previously described. How- 
ever, we also pointed out that the amplitude difference 
is confined to sharp peaks in the time domain. Such a 
confined change in the waveform amplitude leads to a 
Fourier power being dispersed over a large frequency re- 
gion, including outside the LISA band. As a result, there 
is a correspondingly small change in the SNR: of the or- 
der of « 1% relative to the quadrupolar for malism. Such 
a resu lt is in agreement with the analysis of lBabak et al.l 
([20"nl . which was carried out for EMRIs. Therefore, we 
see that the change in SNR is primarily dominated by 
the modifications introduced in the geodesic description 
of the equations of motion, and not in the octopolar cor- 
rection to the waveform generation. 

The analysis presented in this section, in partic- 
ular Figure [5l makes it clear that relativistic cor- 
rections to the waveforms accumulate with multiple 
bursts. In other words, over a single burst (pericen- 
ter passage), a quadrupolar waveform calculation us- 
ing Newtonian dynamics might be sufficient. However, 
if one wants to estimate parameters associated with 
the central MBH, then multiple bursts might be nec- 
essary to associate the events to a single SCO trajec- 
tory. In terms of data analysis, for a detection search 
it is simpler to look for a single burst using tech- 
niques suc h as excess powe r and wavelet decompo sitions 
(e.g. see [Anderson et~aT1 120011: iKlimenko et alJ 12004 
IStuver fc Firm] 120061 : ICamarda fc Ortolanll2006t l For es- 
timating MBH parameters, the results of this paper sug- 
gest that multiple bursts might have to be connected. 
For this to occur, a single template may be used, but as 
our results indicate, the template will need to incorporate 
the effects of general relativity. 10 

At this junction, we should comment on some of the 
caveats in the conclusions derived from the analysis pre- 
sented here. First, in this paper we have primarily con- 
centrated on the question of characterizing the gravita- 
tional waves through the study of the SNR and overlap. 
In this study, however, we have not maximized these data 
analysis measures with respect to (intrinsic or extrinsic) 
orbital parameters. Although it might be possible to im- 
prove the SNR and overlap via parameter maximization, 
our study suggests that the introduction of relativistic 
effects, such as precession, lead to a clear imprint on the 
waveform that might be difficult to mimic with a purely 
Newtonian waveform irrespective of its parameters. Sec- 
ond, in this paper we have only touched the iceberg of the 
signal characterization and parameter estimation prob- 
lem. A possible route to study this problem is through 
a numerical Fisher analysis, with the complications de- 
rived from the fact that the waveforms are known only 
numerically. Furthermore, increasing the complexity of 
the waveform will also increase the computational cost 
of these studies and, thus, it might be interesting to in- 
vestigate whether it is possible or advantageous to search 

10 A note of caution should be added here, since a detailed study 
of the maximization of the overlap with respect to orbital param- 
eters in the Newtonian waveform has not yet been carried out. 
Indeed, it might be possible to mimic some of the relativistic ef- 
fects with Newtonian waveforms with varying parameters, but such 
mimicking is most probably not possible for highly relativistic EM- 
RBs. 



for individual bursts with similar frequency and identify 
them as belonging to the same physical system. These is- 
sues are beyond the scope of this paper, but they should 
be addressed in future investigations. 

Setting these caveats momentarily aside, let us con- 
clude with some discussion of the extreme relativistic 
case introduced earlier. As we have seen, relativistic 
corrections can introduce strong modifications to EMRB 
waveforms, which depend on how relativistic the EMRB 
event is and, in particular, on the pericenter velocity. 
The corrections are particularly strong for the class of 
EMRBs that inhabit the boundary between EMRBs and 
EMRIs, defined by the T cut = 3 x 10 4 s value, corre- 
sponding to the period between apocenter passages. An 
example of such an event is the extreme orbit discussed in 
SectionJU whose waveform is shown in Figure[8l Observe 
that a simplistic Newtonian description misses the rich 
structure, in which the SCO whirls twice about the black 
hole before zooming out to apocenter again. This behav- 
ior is missed entirely when we evolve the orbit with the 
Newtonian equations of motion, even though the same 
initial conditions were used. Even though in the previ- 
ous cases a Newtonian waveform might be able to extract 
relativistic events by adjusting intrinsic parameters, such 
is definitely not the case for the highly-relativistic event 
of Fig.[3 since no choice of parameters in the Newtonian 
waveform can reproduce its multiple-peak structure. 

Due to their whirling behavior, the extreme or- 
bit waveform resembles the zoom-w hirl events oft en 
mentioned in the EMRI litera ture ([Hughes! l2001al fbl: 
iGlampedakis fc Kennefickl [20021. However, the event is 
still an EMRB and not an EMRI because the period be- 
tween apocenter passages is too long. For our galactic 
model, we find the probability of a small region of phase 
space around this orbit to be rather high, 10%. If this 
EMRB is detected with sufficiently high SNR it seems 
plausible that a parameter estimation analysis would al- 
low for a determination of th e background par a meters , 
such as the black hole spin. iBarack fc Cutlen (l2004bf) 
have already investigated LISA's ability to measure MBH 
properties using approximate EMRI signals. They found 
that, depending on the actual orbital parameters, it will 
be possible to measure the MBH spin with fractional er- 
rors of 10~ 3 to 10~ 5 . This high precision measurement 
is the result of observing up to ~ 10 6 complete orbits. 
Conversely, it is very unlikely that EMRB measurements 
will be able to match the measurement capabilities of 
EMRI signals, since only a few bursts will probably be 
available. Whether accurate parameter extraction is pos- 
sible can only be determined with a more detailed data 
analysis investigation of EMRBs. 

6. CONCLUSIONS 

We have studied the effects of relativistic corrections 
on the gravitational waves produced by EMRBs. These 
events originate from long period orbits of a SCO around 
a MBH, leading to large-amplitude, quasi-periodic grav- 
itational wave bursts. Using a more accurate relativis- 
tic treatment of the phenomenon, we have improved on 
the waveforms and trajectories relative to previous work. 
The orbital trajectories were corrected by accounting for 
the spacetime curvature of the system for Schwarzschild 
and Kerr MBHs. The waveform generation was corrected 
by accounting for the next order term in the multipolar 
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Fig. 8. — Plot of the quadrupole Newtonian (solid) and 
Schwarzschild (dashed) gravitational waveform as a function of 
time. 

expansion of gravitational radiation. 

We found that relativistic corrections change the wave- 
form shape relative to its Newtonian counterpart. One 
of the most significant changes was found to be an 
amplitude-modulated dephasing, produced by the rela- 
tivistic corrections to the orbital trajectory and, in par- 
ticular, by relativistic precessional effects. Other effects 
included a change in the amplitude of the waveform, par- 
tially produced by the inclusion of higher-order terms in 
the gravitational wave generation scheme. 

The magnitude of the relativistic corrections was found 
to be directly proportional to the pericenter velocity 
of the orbit, as expected. Surprisingly, we estimated 
that at leas t 50% of the orbits analyzed in iRubbo et al.l 
acquire relativistic velocities and, thus, non- 
negligible relativistic corrections. We investigated these 
corrections in detail by choosing a test orbit, that is rep- 
resentative of the kind tha t dom inated the event rate 
calculation of IRubbo et al.l (|2006l ). We also studied a 
limiting case of a highly relativistic EMRB and found 
that it whirls more than once around the MBH before 
zooming back to apocenter and becoming silent again. 

We have also discussed the possible consequences that 
relativistic effects might have on the detection and pa- 
rameter estimation of gravitational waves from EMRBs 
by LISA, namely a change in the SNR and loss of overlap. 
These changes are mainly due to the relativistic treat- 
ment of the equations of motion, while a quadrupolar 
wave generation formalism seems to suffice. This finding 
is relevant particularly to match filtering searches, where 
a Newtonian treatment of the orbit might lead to a de- 
terioration of confidence limits. Furthermore, our study 
suggests that, given an EMRB gravitational wave detec- 
tion, it might be plausible to extract or bound the spin 



of the central potential with a template that takes into 
account the Kerr character of the MBH. Other astro- 
physical consequences include a possible i ncrease in the 
event rate, w hich implies that the rate s of IRubbo et al.1 
(|2006t 120071 ) and iHopman etlill (|2006h might be lower 
limits, although these estimates are still dominated by 
uncertainties in the astrophysical modeling for the host 
galaxy. 

In addition to the astrophysical modeling, future re- 
search should tackle the details of the EMRB data analy- 
sis and signal extraction issues put forward above. Based 
on the results of this paper, one may explore the possi- 
bility of testing alternative theories of gravity with EM- 
RBs by performing m atched filter i ng with templates from 
alternative theories dWilH Il998t IScharre fc Willi 120021 : 
IWill fc Yun"el2004l ; iBerti et al.ll2005bl |al). Another possi- 
ble avenue for future research is the study of confidence 
limits with which the spin of the central MBH can be 
measured. 

This research could then be used to exam- 
ine whether EMRB events can distinguish between 
a pure Kerr MBH and some other spacctimc. 
Such studies have alr e ady b egan with the analy- 
sis of iCollins fc Hughei pOOik (Glampedakis fc Babald 
(200(1 and lBarausse et al l(|2007D . where comparisons be- 
tween a Kerr and other non-Kerr spacetimes were per- 
formed. Such studies c ould be extended to the pe rturbed 
Kerr solution found by lYunes fc Gonzale3 (|2006f ) , where 
the perturbation is parameterized by the Weyl tensor of 
the external universe and could represent some exter- 
nal accretion disk, planetary system or some other com- 
pact object. Ultimately, these investigations will decide 
whether EMRB events are worth studying in further de- 
tail by future gravitational wave observatories. 
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TABLE 1 

SNR AND OVERLAP FOR DIFFERENT EMRB ORBITS. 



r p [fipc] 


Vp/c 


P [yrs] 


e 


# Bursts 


Ap| 

P 


|A„M| 

pt 1 ) 


(hi\h 2 ) 


(hi\h 2 )W 


7.8 


0.196 


0.0042 


0.603 


6 


0.38 


0.049 


0.27 


0.98 


7.5 


0.207 


U.OUoo 


0.684 


5 


0.13 


O.UOU/ 


0.20 


0.88 


7.2 


0.217 


0.0074 


0.749 


1 


0.22 


0.057 


0.31 


0.90 


7.1 


0.220 


0.0086 


0.777 


3 


0.49 


0.061 


0.33 


0.99 


13 


0.163 


0.0450 


0.859 


1 




0.002 




0.90 


8.0 


0.217 


0.1953 


0.968 


1 




0.033 




0.97 


10 


0.193 


0.7683 


0.984 


1 




0.009 




0.72 


12 


0.173 


3.0407 


0.992 


1 




0.015 




0.80 



Note. — In this table we present the SNR and correlation computed in the frequency 
domain between Kerr quadrupole and octopole waveforms for different EMRB orbits. 
The orbits are separated into two groups: highly-bursting (top) and slowly-bursting 
(bottom). The first five columns present information about the different orbits, including 
how many times they burst, their eccent ricities and periods, which were chosen directly 
from the allowed EMRB phase space of Rubbo ct al. (2006) and thus represent Milky 
Way sources. The sixth and seventh columns present the difference in SNR between a 
Newtonian and Kerr quadrupole waveform relative to the former using the entire data 
set and only one burst respectively. Similarly, the eighth and nineth columns show the 
correlation between the plus polarizations using the entire data set four and one burst 
respectively. Since the slowly-bursting sources burst only once, the sixth and eighth 
columns are redundant for these sources. All calculations assume the observer is located 
on the 2-axis and random initial inclination angles. 
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TABLE 2 

Comparison between quadrupole and octopole waveforms. 



Angle [degrees] a; obs [kpc] y ohs [kpc] z obs [kpc] Amp. diff. 












8 


3.9% 


20 


2.73 





7.51 


3.3% 


40 


5.14 





6.12 


11.5% 


60 


6.93 





4 


15.2% 


80 


7.88 





1.39 


19.7% 


100 


7.88 





-1.39 


21.6% 


120 


6.93 





-4 


17.8% 


140 


5.14 





-6.13 


43.3% 


160 


2.74 





-7.52 


8.6% 



Note. — Here we present an approximate measure of the am- 
plitude difference between the quadrupole and octopole waveforms. 
We concentrate only on the first burst of radiation and we normalize 
the difference to the maximum of the first peak of the quadrupole 
waveform. The difference is presented as a function of the observer 
location, which is always at a fixed radial distance of r D b s = 8 kpc, 
but rotated about the y-axis on x — z plane (6 is here the usual Euler 
polar angle.) 



